Semiclassical quantization of maps with a variable time scale 
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Quantization of energy balance equations, which describe a separatrix - like motion is presented. 
The method is based on an exact canonical transformation of the energy-time pair to the action-angle 
canonical pair, {E,t) {1,8)- Quantum mechanical dynamics can be studied in the framework of 
the new Hamiltonian. This transformation also establishes a relation between a wide class of the 
energy balance equations and dynamical localization of classical diffusion by quantum interference, 
that was studied in the field of quantum chaos. An exact solution for a simple system is presented 
as well. 
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I. INTRODUCTION 

It is known for many applications ||l|,||,H,^ that classical dynamics near hyperbolic points can be described by a 
map, T in the energy-time iE,t) canonical variables T{E^t) —f {E,t). This map is known as the separatrix map, 
and defines a motion in the vicinity of a separatrix, where a period of the unperturbed (twist) map is an arbitrary 
function of energy T{E). Its explicit form is defined by n + 2 turning points as follows 

^ ' \ £;-("-i)/2 if n > 1, ^ ' 

where the energy of the separatrix without loss of generality is taken to be zero and it is assumed that ^ 1. A 
perturbation that is a periodic function of time, for instance, esini^i, that is relevant for a variety of applications, 
may be considered. In this case, T is an energy balance equation that describes the energy change over the period 
T{E): 

En+l = En + eSUWtn 

tn+i^t„+T{En+l), (1.2) 

where v is the frequency of the perturbation of the strength e, and it is assumed that En = E{tn — 0). The period 
of the unperturbed nonlinear motion T{En) describes a wide class of nonlinear systems with variety of applications, 
including: the celestial mechanics of the perturbed Kepler system charge particles in a field of a wave packet 

leading to a separatrix mesh phenomenon in non-KAM (Kolmogorov-Arnold-Moser) systems Q ( withpossible 
applications for atom cooling traps 0), and of electron dynamics in superconducting Josephson junctions Maps 



like (1.2) are also related to a description of Rydberg atoms in a microwave field by the Kepler map |10 111]. Also 
Bloch electr ons in external fields are described by a separatrix map |p^ , p^ , |l5| , p^|jl^ ] and similar systems. 
The map (|1.2|) can be derived from the Hamiltonian, 



H = Ho{E) + V{t, t) = Ho{E) + (e/z/) cos{,^t) ■ 62^^), (1-3) 
where the unperturbed Hamiltonian 
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Ho{E) = ^ y T{E')dE' (1.4) 

depends only on the energy E that is conjugate to the time t, while r is the formal time parameter, and_^2Tr (t) = 
^^_j^(5(t — 27rn) is the periodic ^-function with the period 27r. Quantization of the Hamiltonian (|l.3 ) in the 
framework of the energy-time {E, t) canonical pair has been pres ented in previous studies (see for instance [pyiC JTl| , p^ ). 
The main deficiency of quantization of the Hamiltonian ( |l.3| ) is an appearance of the unphysical time for a wave 
function and this fact has been pointed out earlier, in publications jl^jl^. Classically it is always possible to establish 
a link between the formal time parameter r and the real time t for any individual trajectory by the solution t = t{T) 
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(or by canonical transformation). In quantum mechanics (for a wave function), it is no longer possible. Therefore, 
this approach of quantization in the framework of (-E, t) variables was used to study time independent characteristics, 
such as localization [p| jTc|Jll| ] . On the other hand, properties which are explicitly time dependent, such as dynamical 
correlation functions or dynamical fluctuations, cannot be studied by this approach. Moreover, there is no rigorous 
proof of the validity of these asymptotic solutions. It is also believed that the main deficiency is related to the 
semiclassically approximate definition of the energy operator as E = —ihd/dt, where h is the effective Planck's 
constant. There is no physical justification of this semiclassical definition. Therefore, to overcome these shortcomings 
of the quantum description, it is reasonable to rewrite the system in such a form, where the time parameter appears 
as the physical time. It is convenient to rewrite the system ( |1.2| ) and ( |1.3| ) in terms of the action-angle variables (/, 9), 
related by a canonical transformation H{E, t, r) H(/, 0, t), where t is the real time. This canonical transformation 
must be such that under this variable change one transforms Hq{E) 'Ho{I) and V{t,T) V{0,t). It should be 
stressed that for kicked systems this change of variables has some specific properties, namely, the new potential V{9, t) 
must be independent of / and it must be a function only of the phase 9 and the time t. Otherwise a shift in action 
at a kick is an ambiguous procedure due to the discontinuity of the action I(t — i„) at the n-th kick. 



will be rewritten in the action-angle variables. In 
|16[. Semiclassical quantization of area-preserving 



In what follows, for the semiclassical quantization, the map ( |1.2 
this case the semiclassical quantization procedure is standard |15 

maps was subject of earlier studies. In particular quantization of monotonic twist maps jl^], with interpolating flows 
|l9| was considered. We show first, in Sec. 2, that the map (L2) can be rewritten in terms of action-angle variables. 
Semiclassical quantization and an exactly solvable problem are presented in Sec. 3. The results are summarized in 
Sec. 4. 



II. EXACT CANONICAL TRANSFORMATION {E, T) {I, 0) 



Hamiltonian equations of motion that produce the map (1.2) have the following formal form 

oo 

(a): dE/d9 = -dH/dt = €s\nvt ^ 5{e~2Tm) 



n— — oo 



dt/d9 = dH/dE = T{E) /2Tr, 



(2.1) 



where H is giv en b y ( |l.3| ) and ^ is the formal time-derivative (t = 9). Integration over in the limits (2?™ — 0, 2Trn + 
27r — 0) gives (1.2). We obtain from the equations ( ^.l^ ,b) the Hamiltonian equations for the action-angle variables 
(1,0) and the real time parameter t. Let us start from (p.l|b). Inverting this equation, one obtains 



d0 
'dt 



[T{E)/2n]-'^n{I), 



(2.2) 



where it is supposed that T{E) ^ 0, that means that t is a single valued function of 0. The frequency corresponds 
to some parametric variable change E = E{I) along a trajectory, such that 



dE/dl 



(2.3) 



I n th is case t he e n ergy E — Hoil) is an unperturbed Hamiltonian with the new action variable /, that is defined by 
(2.3). From (2T)-(2_^) we perform the following chain of parametric changes 



dE _dt dE _ dt dE dl 
l£^d0'lt~'(M'~dl''dt 



= I. 



It follows from (p.l|a) and ( |2.4D that 

oo 

i — £ sin i/t S{ 

n— — oo 

where Q{0) is a step function of the following form 



27m) = esini/t— 6(6') , 
o9 



(2.4) 



(2.5) 



11 

27r TT -^-^ I 

1=1 



(2.6) 
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and it is a result of integration of the periodic (5-function written in the form 

oc 11^ 



(2.7) 



1=1 



Equations (^.2[) and (2.5) are the Hamilton equations generated by the Hamihonian 



n = Tioil) - 2e' sin i^t 



CO ^ 



sin/6' 



1=1 



(2.8) 



where e' 



To restore periodicity of the Hamiltonian in 9, the foUowing gauge-hke transformation is carried out 



J = I ^ cos i^t . 

V 



The corresponding Hamiltonian is 

e' °° 1 

n^n^U cos vt) - 2e' sin vtS^ - smW = Ua + V ■ 



(2.9) 



(2.10) 



1=1 



The equations of motion (2.2) and ( |2.5| ) take the form 

J — 2e' sin vi cos IQ 

e = n{j - ^cosvt) . 



(2.11) 



In this form, the equations were derived from the Hamiltonian 7i, therefore the transformation from H to TL, that 
was presented here, is a canonical transformation. 



III. SEMICLASSICAL QUANTIZATION 

A. Floquet operator 

The system (2.11),( pJ^ ) can be simply quantized semiclassically with 

'ihd/dO = hn, 



J ^ J 

where h, in our dimensionless units, is a dimensionless parameter, that plays the role of Planck's constant, while h is a 
quant um numbers operator. In what follows h will be used for the semiclassical quantization. Since the Hamiltonian 
( ^.10 ) is periodic both in time and in the angle 9, the Floquet theory can be used for the analysis. Therefore, a 
solution of the Schrodinger equation 



can be considered due to the Floquet theorem in the following form 

^(0,0 =e-^^*e^'^VA,K(^,i), 
where periodicity of 7i in both 9 and t is taken into account, and the functions 

V'A,« = V'a,k(6', t) = Va,k.(6' + 27r, t + 27r/i/) 



(3.1) 



(3.2) 



(3.3) 



are periodic in time with the periods of the perturbation 2Tr/i' and periodic in 9 with the period 2n, while < k < 1 
and A are a "quasimomentum" and a quasienergy correspondingly. Taking into account that commutation of the 
Hamiltonian with the exponential e*"^ leads to the shift on k for ri. 



nih, 9, i)e'''VA,.(^, t) = e"''>n{n + 9, t)VA,«(e, t), 



(3.4) 
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we can rewrite the Schrodinger equation ( p.l\ ) in the following form 
where F = F{k) is the so-called Floquet operator 

~ d 



(3.5) 



(3.6) 



with A and V'a.k as the eigenvalues and eigenfunctions correspondingly. The solution is cast in the form of the Fourier 
expansion 



(3.7) 



moreover the Fourier harmonics 



\n)\j) 



27r 



(3.8) 



are the eigenfunctions of the unperturbed system with e = 0, and 4>n,j = (i| ('t-|^a,k) arc the coefficients of the 
expansion. Matrix elements of the Floquet operator 



F", ^ {3\{n\F\n')\j' 



(3.9) 



specify the equations for the coefficients (t>n,j- First we calculate the matrix elements for Tig. For this purpose we 
rewrite it formally as 



HoiJ + K - 1 cos vt) = e'v '=°''''*(»a»)7Yo(i + k) , 



(3.10) 



where = d/dn and it is assumed that Tio does not operate on functions of k. Calculation of the matrix elements 
(ph yields for (ETol) 



3,3 
n.n' 



(j1e'^™^''*(*^»)|/)(n|7T!o(J + Ai)|nO =^z™J™(i-aj(5,v,,+,„Ho(n' + K)5„,^ 



(3.11) 



where Jm{x) is the Bessel function. Correspondingly, the matrix elements for the additive part of the perturbation 
in ( |2lc| ) are 



y3:3 
n.n 



\^\-^3',3 + l-^3',3-l\ ■ "7^*^' 



(3.12) 



These matrix elements of ( 3.11 ) and ( 3.12| ) lead to the following equation for the expansion coefficients 
e' e' 1 

V'r"J„(i — (9K)Ho(n + K)cj)n,j+ni + 77 ri'/'n'.j + l " 4>n\j-l\ = + ^i)(t>n,j ■ (3.13) 

m n'^n 

It is similar to a quasi-lD Anderson-like chain for the cases when ^{^{n + k) corresponds to a random potential. It 
i s al so the generalization of a specific application in pO| . The classical motion is chaotic in many types of T{E) in 
(1.2). Therefore, the obtained equation establishes a relation between a wide class of energy balance equations of the 
form (1.2), (1.3) and dynamical localization of the classical chaotic diffusion In the general case this equation is 
useful for numerical studies. It is analytically tractable in perturbation theory for the parameter ^ <C 1. Therefore, 
it is instructive to present an exactly solvable example where (1.4) is linear, with a constant period T{E) — 27r/aj. In 
this case a solution can be obtained in an explicit form analytically. 
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B. An exact model 



An exact solution can be obtained for the harmonic osciUator Ti.o{I) — to I with constant frequency = lo. The 
hnear driven system described by the Hamihonian Ti{I,6,t) is integrable. The equation ( |3.13 ) takes the foUowing 
simple form 



For the specific case of the harmonic oscillator k = 0. For other linear models, for example models of an appropriate 
band structure in solids, k, does not have to vanish. A solution of equation (3.14) is cast in the form of a sine function 
by the following substitution 



4>n.j = .^jSinc [^{n — X/lo + k)]. 



(3.15) 



where the sine function is sine x = and Zj has to be determined. In this ease the quasienergy spectrum turns 
out to be 



A = LOK + Lon . 



Inserting (B.6) in (3.14) and taking into account that 

sinc7r(n' — n) = 0, 



(3.16) 
(3.17) 



(3.18) 



since sin7r(n' — n) = 0, one obtains from ( 3.14 ) the following relation for the Zj 

2jZj + z[Zj+i + = 0, 

where z = e'u/hv'^. The solution is (22) 

Z,^Z,{z) = {-iyj,{z). (3.19) 

The result ( 3.16| ) and ( 3.19 ) is only a specific solution. The general solution is found by replacing in ( |3.14 ) j by 
i — jo and A by A — vj^ . The resulting solution of the eigenvalue problem is 

Xjg = i/jo + UK + urn (3.20) 

and the corresponding eigenfunction is 

Z, = i-iy-'"Jj-j,iz). (3.21) 



Floquet theory implies that only Xj^ in an interval of size i^, say [0, i^) should be used (see (3.2)). In particular for the 
harmonic oscillator 



(3.22) 



if it is in the interval [0, v), as expected. 



IV. SUMMARY 



Quantization of energy balance equations, in the form of the map (1^), is presented. This procedure of quantization 
consists of the two steps. The first one is a transformation of the Hamiltonian (1.2) to the form (2.10) by the canonical 
transformation {E,t) — )■ {1,0). The second stage is semielassical quantization. The Hamiltonian formulation of the 
problem allows semielassical quantization of action-angle variables in the framework of the standard procedure: 
1—^1 = —ihd/dO = hh, where ft, is a dimensionless semielassical parameter such that the classical limit is ft. ^ 0. 
It should be noted that the semielassical approximation requires that the width of the potential V is larger than 
the de-Broglie wavelength. Therefore, it is necessary to truncate the Fourier expansion in (^^, corresponding to 
replacement of the (5-function by a function of finite width. The results obtained in Sec. 3.2 in the semielassical 
approximation are correct only up to the leading correction in h to the classical ones [|l5|. The advantage of the 
present analysis is that the transition from (1.2) to (2TC) is exact, and wave functions depend on physical time, while 
the main deficiency of straightforward quantization of (1.2) is an appearance of an unphysical time parameter |p^ , pl| . 
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